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This  Report  describes  a  theoretical  method,  based  on  linear  stability 
theory,  for  predicting  flow  transition  on  axisymmetric  heated  underwater 
bodies,  and  a  simple  experiment  using  a  four-to-^ne  ellipsoidal  nose  in  the 
30-inch  Water  Tunnel. at-A5ffi~fTeddingtt!Tr)  ‘  The  experiment  showed  that 
filtering  the  tunnel  water  and  vibration  of  the  model  have  an  effect  on 
transition. 

Agreement  between  theory  and  experiment  is  good  without  heating,  but  a  more 
sophisticated  experiment  is  needed  to  validate  the  theory  for  heated 
bodies.  Although  increases  in  transition  Reynolds  number  of  40  per  cent 
were  observed  experimentally,  the  theory  suggests  that  much  greater 
increases  are  possible  for  the  same  heat  input  by  distributing  it 
differently  over  the  body  surface,  lu  .-  .  C ■  s  *  r  ■  r 


(An  earlier  version  of  this  Report  was  presented  orally  at  Euromech  228 
Colloquium  on  Boundary  Layer  Instability  and  Transition,  Exeter,  September 
1987). 
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Notation 


D  Body  diameter  (m) 

Cp  Specific  heat  of  water  (J  kg'1  'C'1) 
g  Acceleration  due  to  gravity  (ms"2) 

T  Temperature  (‘C) 

u  Fluid  velocity  (ms'1) 

x  Axial  distance  from  nose  (m) 

B  Coefficient  of  thermal  expansion  of  water  (‘C  *) 
v  Kinematic  viscosity  (m2s"1) 

Subscripts 

w  Denotes  value  at  body  surface 
»  Denotes  ambient  or  freestream  conditions 


v 


TRANSITION  CONTROL  VIA  BOUNDARY  LAYER  HEATING 


By  D  J  Atkins 
S  J  Pearce 

1.  INTRODUCTION 

It  is  well-known  that  the  transition  Reynolds  number  on  flat  plates,  and 
axisymmetric  bodies  with  a  favourable  pressure  gradient  along  most  of  their 
length,  can  be  significantly  increased  by  heating  the  surface.  This 
phenomenon  has  implications  for  drag  reduction  and  has  been  the  subject  of 
much  published  research,  particularly  in  the  USA.  Wazzan  et  al 
(References  1,  2)  showed  theoretically,  for  a  flat  plate  heated  to  a 
constant  temperature  above  ambient,  that  the  critical  Reynolds  number 
(based  on  displacement  thickness)  for  neutral  stability  was  increased  from 
520  to  approximately  15000.  The  optimum  overheat  was  about  40°C  giving  a 
transition  Reynolds  number  of  about  2  x  108  as  compared  with  the  unheated 
value  of  3  x  106.  These  findings  were  confirmed  for  small  overheats  up  to 
8°C  by  Strazisar  et  al  (Reference  3).  Barker  and  Jennings  (Reference  4) 
and  Barker  and  Gile  (Reference  5)  investigated  the  stability  of  a  heated 
laminar  boundary  layer  developing  inside  a  cylindrical  tube.  A  number  of 
different  wall  temperature  distributions  were  tried  and  the  optimum 
increase  of  a  factor  of  3-4  in  transition  Reynolds  number  was  observed  for 
a  small  (7-8°C)  uniform  overheat.  Additional  heating  showed  no  further 
increases,  in  contradiction  to  the  theory,  which  predicted  increases  up  to 
an  overheat  of  35°C.  Lauchle  and  Gurney  (Reference  6)  studied  a  "laminar- 
flow"  body  with  a  large  region  of  favourable  pressure  gradient  which  was 
heated  electrically.  The  maximum  heating  power  available  was  93  kw,  which 
increased  the  transition  Reynolds  number  from  4-5  x  10°  to  3-6  x  107.  The 
effect  of  increasing  the  total  heating  power  from  zero  to  the  maximum  value 
resulted  in  a  linear  increase  in  transition  Reynolds  number  at  first, 
followed  by  a  levelling  out  at  higher  values.  The  theory  however  predicted 
a  linear  increase  throughout  the  entire  range  of  power  inputs.  The 
disparity  between  theory  and  experiment  was  attributed  by  the  authors  to 
surface  roughness  and  the  effects  of  free  stream  particles,  and  further 
research  has  been  carried  out  on  these  effects  (References  7,  8).  This 
showed  that  heated  boundary  layers  are  more  susceptible  to  external 
disturbances  than  unheated  ones  and  that  any  improvements  arising  from 
surface  heating  can  be  nullified  by  such  effects. 

All  the  workers  mentioned  above  used  a  theoretical  approach  based  on  two- 
dimensional  linear  stability  theory.  This  involves  solving  a  modified 
version  of  the  well-known  Orr-Sommerfeld  equation,  in  which  only  the 
variation  of  viscosity  (as  a  function  of  temperature)  is  taken  into 
account.  Lowell  and  Reshotko  (Reference  9)  derived  the  full  disturbance 
equations  in  which  fluctuations  in  all  the  fluid  properties  are  taken  into 
account,  giving  a  more  complex  sixth-order  system,  and  they  showed  that 
there  is  no  major  quantitative  difference  between  the  results  obtained 
using  the  fourth-  and  sixth-order  systems.  However  they  neglected  buoyancy 
effects  in  their  linear  stability  model,  and  in  common  with  previous 
workers  assumed  that  buoyancy  effects  in  the  mean  flow  are  negligible. 
Barker  and  Gile  (Reference  5)  and  Lauchle  &  Gurr»y  (Reference  6)  commented, 
after  some  discussion,  that  the  disagreement  between  their  experiments  and 
the  theory  cannot  be  attributed  to  flow  asymmetries  arising  from  buoyancy 
effects.  Lauchle  and  Gurney  had  three  hot-film  probes  at  the  same  axial 
location,  but  at  120*  intervals  round  the  circumference  of  their  test  body, 
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to  measure  any  circumferential  variations  in  stability.  However  two  of  the 
probes  failed  to  operate  and  they  were  only  able  to  use  the  results  from 
the  probe  on  top  of  the  body.  Asrar  (Reference  10)  showed  theoretically  by 
doing  a  full  three-dimensional  linear  stability  analysis  (again  neglecting 
buoyancy  effects)  that  two-dimensional  disturbances  are  the  most  critical 
for  heated  flat  plate  and  axisymmetric  flows.  He  also  investigated  the 
effects  of  different  surface  temperature  distributions  on  a  6  to  1  prolate 
spheroid  and  concluded  that  a  gradually  increasing  overheat  with  downstream 
position  was  better  than  a  uniform  or  gradually  reducing  overheat. 

The  objectives  of  the  present  research  are  to  develop  a  theoretical 
transition  prediction  method  for  underwater  axisymmetric  heated  bodies  and 
validate  it  experimentally.  The  prediction  method  was  developed  partly  in- 
house  and  partly  under  contract  with  BMT  Ltd  under  the  direction  of 
Dr  M  Gaster.  Two  computer  programs  known  as  HEATEDBL  and  FLTRANS  were 
written  and  these  are  described  in  Reference  11.  Brief  details  of  the 
experiment  are  given  in  Section  3  and  more  detailed  descriptions  are  given 
in  Appendix  A  and  Reference  12. 

2.  THEORY 

The  prediction  method  involves  computing  the  laminar  boundary  layer  and 
linear  spatial  stability.  The  boundary- layer  calculation  solves  the 
coupled  momentum  and  energy  equations  on  a  grid  which  expands  with  the 
boundary  layer  based  on  the  usual  similarity  transform  approach.  The 
fourth-order  Orr-Sommerfeld  equation  with  variable  viscosity  is  solved 
using  a  formulation  of  the  compound  matrix  method  developed  at  BMT  Limited 
by  Willis  (Reference  13).  Following  Smith  and  Gamberoni  (Reference  14), 
the  onset  of  transition  is  assumed  to  occur  where  the  most  unstable 
frequency  is  amplified  by  a  factor  en,  where  n  usually  takes  the  value  9. 
Buoyancy  effects  have  been  neglected;  a  discussion  in  Appendix  B  suggests 
that  for  the  present  geometry  and  flow  conditions  they  are  not  significant. 

Numerical  experiments  on  flat  plates  carried  out  by  Atkins  (Reference  15) 
have  shown  that  the  compound  matrix  method  requires  many  data  points 
(typically  between  400  and  3200  depending  on  the  overheat)  to  produce 
accurate  eigenvalues,  but  is  relatively  insensitive  to  other  physical  and 
numerical  parameters.  Calculations  for  the  temporal  stability  case  using 
the  full  sixth-order  equation,  which  takes  buoyancy  effects  into  account  as 
well  as  fluctuations  in  all  the  other  fluid  properties,  have  been  carried 
out  by  Thomas  (Reference  16).  These  have  shown  non-trivial  but  small 
differences  between  fourth-order  and  sixth-order  eigenvalues,  confirming 
the  results  of  Lowell  and  Reshotko  (Reference  9). 

3.  EXPERIMENTAL  INVESTIGATION 

The  30-inch  Water  Tunnel  at  ARE  (Teddington)  was  used  for  a  simple  heated 
body  experiment.  The  test  body  was  a  four-to-one  ellipsoidal  nose  attached 
to  a  long  cylinder.  Figure  1  shows  a  sketch  of  the  model  geometry  and  the 
pressure  coefficient  which  was  calculated  using  the  method  of  Hess  and 
Smith  (Reference  17).  Body  offsets  and  pressure  coefficients  are  tabulated 
in  Reference  11.  Unlike  previous  body  shapes  tested  experimentally,  the 
current  body  has  a  region  of  mildly  adverse  pressure  gradient  which  is 
characteristic  of  typical  underwater  vehicle  shapes.  Heating  was  supplied 
by  pumping  hot  water  through  a  copper  coil  wound  helically  inside  the  nose. 
Intermittency  measurements  were  made  at  a  number  of  axial  positions  using 
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glue-on  hot-film  probes.  Body  surface  temperature  was  measured  at  the  aft 
end  of  the  heating  coil  using  a  thermocouple  and  at  a  few  other  locations 
by  calibrating  the  hot-film  probes  as  thermocouples.  The  temperature  data 
from  the  hot-film  probes  turned  out  to  be  rather  unreliable.  The  model 
surface  was  made  smooth  enough  to  avoid  premature  tripping  due  to 
roughness,  and  buoyancy  effects  were  thought  to  be  insignificant.  However, 
the  effects  of  model  vibration  and  filtering  the  tunnel  water  were  found  to 
be  significant. 

The  turbulence  intensity  was  measured  in  the  empty  tunnel  and  the  results 
are  shown  in  Figure  2.  There  is  a  honeycomb  arrangement  in  the  settling 
section  of  the  tunnel  upstream  of  the  contraction,  and  this  results  in  low 
turbulence  levels  of  0-1  per  cent  at  the  start  of  the  working  section  over 
a  wide  range  of  tunnel  speeds  (3-13  ms'1).  However  the  turbulence  levels 
increase  significantly  further  downstream,  owing  to  the  effect  of  the 
slotted  wall  section.  As  a  result  of  this,  the  model  was  placed  as  far 
forward  as  possible,  with  the  nose  being  0-69  m  from  the  tunnel  throat. 

The  range  of  tunnel  speeds  obtained  during  the  transition  experiments 
ranged  from  2  to  14  ms'1.  At  speeds  above  10  ms'1  there  was  evidence  of 
low-frequency  model  vibration  which  was  clearly  visible  at  14  ms'1. 

Typical  signals  from  the  hot-film  gauges  are  shown  as  Figure  3.  As  in 
previous  work  (Reference  18),  four  different  flow  regimes  are  identifiable, 
viz:  laminar,  wave-like,  intermittent  and  fully  turbulent.  Figures  4(a) 
and  4(b)  show  the  spectral  density  of  the  wave-like  signals.  The  unheated 
spectrum  has  a  single  peak,  which  suggests  that  the  disturbance  may  be  a 
Tollmein-Schlichting  wave  (Reference  19).  The  heated  spectrum  is  more 
complicated  with  a  fundamental  peak  and  first  and  second  harmonics,  which 
suggests  that  there  may  be  non-linear  effects. 

Two  separate  series  of  transition  experiments  were  carried  out.  The  first 
series  was  without  heating  and  the  second  mainly  with  heating.  Halfway 
through  the  second  series,  the  tunnel  water  was  filtered  for  eight  hours 
and  this  had  an  effect  on  the  results  (see  Section  5).  With  heating, 
increases  in  transition  Reynolds  numbers  of  at  least  40  per  cent  were 
observed,  and  intermittency  values  were  generally  lower  than  the 
corresponding  unheated  values.  However  the  heated  boundary  layer  was 
observed  to  be  more  susceptible  to  external  effects  such  as  model 
vibration,  the  presence  of  particles  in  the  flow,  and  air  or  vapour 
bubbles. 

Because  of  the  unreliability  of  the  hot-film  temperature  measurements,  the 
variation  of  surface  temperature  (or  surface  heat  flux)  over  the  heated 
portion  of  the  body  surface  was  to  a  large  extent  unknown.  The 
measurements  (Reference  12)  suggest  that  the  temperature  rises  to  a  maximum 
and  falls  off  near  the  tail  end.  The  discussion  in  Appendix  B  suggests 
that  three-dimensional  effects  are  negligible  on  the  present  geometry.  The 
total  heating  power  supplied  to  the  model  was  estimated  from  the  water 
temperature  at  the  coil  inlet  and  outlet  and  found  to  vary  between  4  and 
12  kW  over  the  range  of  flow  speeds  (4-14  ms'1)  and  reservoir  tank 
temperatures  (48-68°C  or  31-51 "C  above  ambient). 

4.  NUMERICAL  CALCULATIONS 

Numerical  calculations  were  carried  out  on  the  unheated  body  for  a  range  of 
flow  speeds  from  4  to  14  ms'1  using  computer  programs  HEATEDBL  and  FLTRANS 
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(Reference  11).  Calculations  were  also  carried  out  on  the  heated  body 
(assuming  axial  symmetry)  at  flow  speeds  of  8,  10  and  12  ms-1  with  a  view 
to  gaining  some  knowledge  of  the  influence  of  surface  heat-flux  variation 
on  transition  to  turbulence  for  the  ellipsoidal  nose.  Heat  was  applied  to 
the  same  region  of  the  surface  as  used  in  the  experiment.  Various  simple 
heat-flux  profiles  were  tried,  viz:  constant,  linearly  increasing  and 
linearly  decreasing.  A  few  computations  were  carried  out  using  profiles  of 
heat  flux  proportional  to  ±Jx  and  ±x2,  rather  than  ±x,  but  the  results 
obtained  were  not  significantly  different. 

Figure  5  shows  surface  temperature  profiles  for  a  total  heat  input  of 
1-25  kW.  In  calculating  these  profiles,  the  effects  of  heat  conduction 
along  the  surface  were  not  taken  into  account.  The  temperature 
measurements  estimated  by  calibrating  the  hot-film  probes  as  thermocouples 
indicated  that  the  temperature  decreased  downstream  of  x/D=20.  This 
suggests  that  the  theoretical  case  of  linearly  decreasing  heat  flux 
corresponds  most  closely  to  the  experiment.  Figure  6  shows  the  values  of 
excess  temperature  at  the  end  of  the  heating  coil,  i.e.  at  x/D=2-2,  plotted 
against  the  total  heating  power. 

Figure  7  shows  predictions  of  flow  transition  based  on  e9  obtained  at  a 
flow  speed  of  10  ms'1  with  constant,  linearly  increasing  and  linearly 
decreasing  heat-flux  profiles.  If  the  total  heating  power  is  less  than 
0-5  kW,  there  is  little  difference  between  the  results  using  different 
profiles.  Above  0-75  kW  however,  different  heat-flux  profiles  give 
significant  differences.  The  results  indicate  that  gradually  increasing 
heat  flux  has  a  more  beneficial  effect  on  transition  than  a  constant  heat 
flux  which,  in  turn,  has  a  more  beneficial  effect  than  a  gradually 
decreasing  heat  flux.  To  illustrate  the  effect  of  flow  speed,  the  results 
for  8  and  12  ms'1  are  shown  for  the  case  of  linearly  decreasing  heat  flux. 
The  trends  in  the  present  results  are  similar  to  those  found  by  Asrar 
(Reference  10). 

5.  COMPARISON  OF  THEORY  WITH  EXPERIMENT, 

5.1.  Unheated  body 

The  unheated  body  experiment  has  resulted  in  a  set  of  consistent  repeatable 
data  which  can  be  used  for  comparison  with  the  theory.  The  measured  flow 
regimes  are  compared  with  theoretical  disturbance  amplification  factors  in 
Figure  8.  Figure  9  shows  the  measured  and  predicted  spectral  density  peak 
frequencies  which  correlate  well.  Figure  10  shows  the  variation  of 
intermittency  with  axial  position.  The  uncertainty  bands  indicate  the 
spread  of  the  experimental  data  over  different  runs.  The  results  at  a 
tunnel  speed  of  14  ms’1  are  somewhat  inconsistent  with  the  rest  of  the 
data,  but  the  model  was  affected  significantly  by  vibration  at  this  speed. 
The  effect  of  filtering  the  tunnel  water,  shown  in  Figure  11,  results  in  a 
small  increase  in  transition  Reynolds  number.  Figures  10  and  11  indicate 
that  transition  is  not  instantaneous  but  occurs  gradually  over  a  finite 
longitudinal  distance.  The  extent  of  the  transition  region  was  estimated 
crudely  by  calculating  the  approximate  slope  of  the  intermittency  curve  at 
an  intermediate  value  and  extrapolating  to  0  per  cent  and  100  per  cent. 
Figure  12  shows  the  extent  of  the  transition  region  estimated  in  this  way, 
which  is  generally  in  the  range  0-25D  -  0-40D  (0055m  -  0-085m).  The  onset 
of  transition,  corresponding  to  the  start  of  turbulent  bursting,  appears  to 
correspond  with  e9  at  low  tunnel  speeds  and  e7  at  high  speeds.  Pearce  and 
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Turner  (Reference  12)  suggested  that,  at  high  tunnel  speeds,  transition  is 
affected  by  model  vibration.  Power  (Reference  18)  has  conducted  an 
experiment  over  a  similar  Reynolds  number  range  on  a  3-6  to  1  ellipsoidal 
nose  and  observed  that  the  onset  of  transition  corresponded  to  between  e8 
and  e12  depending  on  Reynolds  number.  For  other  body  shapes  over  the  same 
Reynolds  number  range,  he  obtained  values  between  e7  and  e12.  Power's 
experiments  were  conducted  in  a  towing  tank  which  probably  had  a  lower 
level  of  turbulence  than  the  ARE  Water  Tunnel. 


In  boundary-layer  and  drag  calculation  methods  which  assume  that  transition 
is  instantaneous  and  require  a  value  for  transition  location  to  be 
specified,  the  experimental  value  corresponding  to  50  per  cent 
intermittency  is  often  used.  These  values  are  indicated  on  Figure  12  and 
the  effect  of  filtering  the  tunnel  water  is  also  shown.  The  transition 
Reynolds  numbers  Us/i/,  based  on  tunnel  speed  and  surface  arc  length,  were 
calculated  at  50  per  cent  intermittency.  The  average  increase  as  a  result 
of  filtering  was  about  2  per  cent. 

5.2.  Heated  body 

Because  of  the  lack  of  reliable  information  regarding  the  temperature  or 
power  variation  over  the  heated  region,  no  direct  comparison  of  theory  and 
experiment  is  possible.  However,  an  attempt  to  correlate  theory  and 
experiment  is  desirable  with  a  view  to  making  recommendations  for  a  future 
experiment. 

In  the  experiment,  the  total  power  supplied  to  the  heating  coil  was 
calculated,  assuming  a  flow  rate  of  2  gallons  per  minute,  according  to  the 
formula 


Heating  power  (Watts)  =  0-15  Cp  (Tin-Tout) 


where 


Cp  =  4190  J  kg-1  °C_1 

Figure  13  shows  the  total  power  supplied  to  the  heating  coil  at  tunnel 
speeds  of  8,  10  and  12  ms'1,  which  is  plotted  against  the  excess 
temperature,  T  -  T(ambient),  as  measured  by  the  thermocouple  embedded  in 
the  model  skin  at  x/D  =  2-2.  There  is  considerable  scatter  in  the  data  and 
so  straight-line  least  squares  fits  were  obtained  which  are  the  broken 
lines.  The  continuous  lines  are  theoretical  curves  of  the  total  power  used 
to  heat  the  surface  over  the  coil  region  plotted  against  the  excess 
temperature  on  the  surface  at  the  same  location.  A  linearly  decreasing 
heat  flux  profile  was  assumed.  In  order  to  obtain  the  same  temperature  as 
measured  experimentally  at  the  thermocouple  location,  the  theoretical  power 
used  to  heat  the  surface  is  only  a  small  fraction  (calculated  to  be  between 
4  per  cent  and  10  per  cent)  of  the  total  power  supplied  to  the  coil  in  the 
experiment. 

Figure  14  shows  theoretical  curves  of  the  variation  of  axial  position, 
corresponding  to  an  amplification  factor  of  e9,  with  power  used  to  heat  the 
surface  at  flow  speeds  of  8,  10  and  12  ms'1.  Once  again,  a  linearly 
decreasing  heat  flux  profile  was  used.  These  results  indicate  that  a  power 
input  of  1  kW  should  give  an  increase  in  transition  Reynolds  number  of 
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about  40  per  cent  at  all  three  flow  speeds.  The  relatively  few  heatei 
measurements  available  are  listed  in  Table  1  and  these  show  increases  in 
transition  Reynolds  number  of  between  17  per  cent  and  40  per  cent.  The 
total  power  supplied  to  the  heating  coil  varied  from  4-7  kW  to  10-9  kW. 

The  measurements  of  intermittency  (Reference  12)  suggest  that  filtering  the 
tunnel  water  has  a  greater  effect  on  transition  on  a  heated  body,  as 
compared  with  an  unheated  body.  A  direct  comparison  between  unfiltered  and 
filtered  water  was  only  possible  at  one  condition  (10  ms*1  and  5  kW  input 
power)  and  this  gave  an  increase  of  about  6  per  cent  in  transition  Reynolds 
number. 

Figure  14  also  shows  the  unheated  body  measured  transition  locations  (the 
circles  indicating  the  experimental  points  correspond  to  50  per  cent 
intermittency).  The  theoretical  power  inputs  which  would  give  the  same 
measured  increase  in  transition  Reynolds  number  are  also  shown  in  Table  1 
and  Figure  14.  These  values  range  from  7-13  per  cent  of  the  total  power 
supplied  to  the  heating  coil.  It  was  possible  to  estimate  the  peak 
frequency  of  the  Tollmein-Schlichting  waves  for  some  of  the  measurements, 
and  these  are  compared  with  theoretical  predictions  in  Figure  15.  It  was 
assumed  that  10  per  cent  of  the  power  supplied  to  the  coil  was  used  to  heat 
the  surface.  The  measurements  and  predictions  are  reasonably  consistent, 
although,  at  the  higher  power  input,  the  measured  peak  frequencies  are 
higher  than  predicted. 

The  previous  discussion  suggests  that  the  theoretical  power  inputs  required 
to  achieve  the  measured  thermocouple  temperatures  are  similar  to  those 
giving  the  measured  increases  in  transition  Reynolds  number.  This  implies 
that  in  the  experiment,  only  a  small  amount  (5-10  per  cent)  of  the  heat 
supplied  to  the  coil  actually  reached  the  portion  of  the  body  surface  where 
the  coil  was  located.  One  possible  explanation  for  this  is  that  there  may 
have  been  significant  heat  conduction  along  the  model  surface  towards  the 
nose  or  tail  regions.  Some  numerical  calculations  were  carried  out  with 
heat  supplied  to  the  nose  region  as  well  as  the  coil  region,  and  these 
showed  that  the  extra  heat  supplied  to  the  nose  region  has  a  relatively 
small  effect  on  transition,  because  the  unheated  boundary  layer  is  stable 
in  this  region.  In  view  of  the  simplicity  of  the  experiment  and  the 
limited  instrumentation,  further  discussion  of  the  results  is  of  limited 
value.  A  detailed  quantitative  validation  of  the  theory  clearly  requires  a 
much  more  sophisticated  and  better  instrumented  model. 

6.  CONCLUSIONS  AND  FURTHER  RESEARCH 

A  theoretical  transition  prediction  capability,  based  on  linear  stability 
theory,  has  been  developed  for  axisymmetric  underwater  bodies,  which 
includes  the  effects  of  surface  heating.  The  30-inch  Water  Tunnel  at 
ARE  (Teddington)  has  been  used  for  a  simple  transition  experiment  involving 
measurements  of  intermittency  using  hot-film  gauges.  The  experiment  showed 
that  transition  was  affected  by  model  vibration  and  filtering  the  tunnel 
water  (more  so  with  heating  than  without  heating).  Heating  the  body 
surface  also  appeared  to  modify  the  Tollmein-Schlichting  wave  spectrum  in 
the  way  suggested  by  the  theory. 

In  the  absence  of  heating,  the  measured  transition  locations  corresponding 
to  the  start  of  turbulent  bursting  agree  with  predictions  of  peak 
amplification  ratios  of  e®  at  low  tunnel  speeds  and  e7  at  high  speeds.  The 
difference  at  high  speeds  is  thought  to  be  due  to  model  vibration.  The 
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signals  from  the  hot-film  gauges  indicate  the  existence  of  Tollmein- 
Schlichting  waves  and  the  predicted  and  measured  spectral  peak  frequencies 
agree  well. 

Numerical  calculations  for  the  heated  body  using  different  heat  flux 
profiles  have  shown  that  a  linearly  increasing  profile  has  a  more  favour¬ 
able  effect  on  transition  than  a  corstant  profile  which  is,  in  turn,  more 
favourable  than  a  linearly  decreasing  profile.  The  measurements  suggest 
that  the  latter  profile  most  closely  matches  the  experimental  situation. 
Comparison  of  theory  and  experiment  suggests  that  only  a  small  fraction 
(about  10  per  cent)  of  the  power  supplied  to  the  coil  was  used  to  heat  the 
surface.  The  theory  suggests  that  much  greater  increases  in  transition 
Reynolds  number  are  possible  than  those  measured  experimentally 
(17-40  per  cent)  by  improving  the  efficiency  of  the  heating  system  and 
distributing  the  heat  differently  over  the  model  surface. 

The  30-inch  Water  Tunnel  can  be  used  for  transition  experiments,  but  it  is 
desirable  to  reduce  the  turbulence  levels  in  the  working  section.  A  new 
closed  working  section  can  be  fitted  which  achieves  this.  The  theory  has 
beer,  validated  quantitatively  for  unheated  bodies  and  qualitatively  for 
heated  bodies.  In  order  to  achieve  a  quantitative  validation  of  the  theory 
for  heated  bodies,  a  more  sophisticated  experiment  with  controllable  heat 
distribution  and  detailed  and  accurate  temperature  measurement  over  the 
whole  body  surface  is  required.  Further  investigations  of  the  receptivity 
of  heated  boundary  layers  to  external  disturbances  and  of  non-linear 
effects  are  also  needed. 
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HEATED  BODY  TRANSITION  MEASUREMENTS 
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Appendix  A 


DESCRIPTION  OF  THE  EXPERIMENT 


A  survey  of  the  available  Test  Facilities  indicated  that  the  30-inch  Water 
Tunnel  at  ARE  (Teddington)  (Reference  20)  would  be  the  most  suitable  for  a 
preliminary  transition  experiment,  because  of  the  relatively  high  Reynolds 
numbers  obtainable.  The  tunnel  has  a  slotted  wall  working  section  4-58  m 
long  and  0-76  m  in  diameter  designed  to  accommodate  relatively  large  models 
(up  to  0-25  m  in  diameter),  and  flow  speeds  of  up  to  20  ms'1  can  be 
achieved.  The  objectives  of  the  experiment  were  to  determine  whether  or 
not  the  30- inch  Water  Tunnel  is  suitable  for  transition  experiments,  to 
develop  suitable  instrumentation  techniques  for  measuring  transition  on 
heated  bodies,  to  determine  the  effect  of  small  increases  in  surface 
temperature  on  transition  for  a  body  with  a  mildly  adverse  pressure 
gradient,  and  to  provide  data  for  an  initial  validation  of  the  theory. 

The  four-to-one  ellipsoidal  nose  (see  Figure  1)  was  specially  manufactured 
and  connected  to  an  existing  cylindrical  afterbody.  The  joint,  at  0-518  m 
from  the  nose  (x/D  -  2-4),  was  hydrodynamically  rough,  which  meant  that  no 
measurements  could  be  carried  out  downstream  of  this  position.  Profile 
tolerances  were  ±  0-0001  m  with  a  limit  of  waviness  slope  of  0-001. 
Hendricks  and  Ladd  (Reference  7)  have  shown  that  the  surface  finish  is  more 
critical  for  a  heated  body  experiment  than  for  an  unheated  experiment,  and, 
based  on  their  results,  a  surface  finish  of  0-2  p.m  rms  was  used. 

The  nose  was  heated  by  pumping  hot  water  around  a  copper  coil  of  0-01  a 
internal  diameter  wound  helically  inside  the  model.  The  extent  of  the  coil 
was  from  0-162  to  0-476  m  from  the  nose  (x/D  -  0-75  to  2-20).  The  gap 
between  the  coil  and  the  inner  surface  of  the  model  was  filled  with  low 
melting  point  Serrabend  (a  lead/aluminium  alloy)  to  improve  thermal  contact 
between  the  surfaces.  Hot  water  was  supplied  from  an  external  tank  fitted 
with  three  3  kW  immersion  heaters  and  pumped  from  the  front  to  the  back  of 
the  model  at  a  nominal  flow  rate  of  2  gallons/min  (0-15  kg  s'*).  For  one 
run,  the  direction  of  flow  was  reversed.  Intermittency  was  measured  using 
DANTEC  type  55R47  glue-on  hot-film  gauges  positioned  at  x/D  -  1-5,  1-75, 
20,  2-2  and  2-35.  These  were  used  rather  than  the  more  sophisticated 
flush-mounted  probes  to  simplify  the  model  manufacturing  process.  Because 
the  gauges  were  designed  primarily  for  use  in  air,  they  were  given  added 
protection  by  coating  with  an  air-drying  polyeurethane  coating  and  the 
leading  edges  of  each  foil  were  trimmed  to  avoid  premature  tripping  of  the 
boundary  layer  by  the  gauge.  Temperature  was  measured  using  K-type 
(NiCr/NiAl)  thermocouples  which  were  positioned  in  the  tank,  the  coil  inlet 
and  outlet  and  the  model  shell  0-001m  from  the  outer  surface  at  x/D  -  2-2. 
The  ambient  temperature  of  the  water  in  the  tunnel  was  also  measured  and 
found  to  be  relatively  constant  having  a  value  of  approximately  17”C.  The 
hot- film  gauges  at  x/D  -  2-0,  2-2  and  2-35  were  also  calibrated  as 
thermocouples,  but  the  repeatability  of  the  temperature  measurements  so 
obtained  was  poor. 
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Appendix  B 


ESTIMATION  OF  BUOYANCY  EFFECTS 


It  Is  Important  to  consider  buoyancy  effects  In  the  laminar  boundary  layer 
on  a  heated  underwater  body.  Although  no  study  appears  to  have  been  made 
on  conventional  vehicle  shapes ,  Yao  et  al  (Reference  21)  studied 
theoretically  the  development  of  a  three-dimensional  water  boundary  layer 
along  a  heated  longitudinal  cylinder.  They  concluded  that  near  the  leading 
edge,  the  boundary - laye r  flow  is  forced- convection  dominant,  but  further 
downstream,  it  becomes  free -convection  dominant.  Their  numerical  results 
for  uniform  surface  temperature  show  that  buoyancy  effects  are  negligible 
when 


x  <  0-1  a/7« 
and  small  up  to 
x  -  0-3  a/7e 


(Al) 


(A2) 


where  a  is  the  cylinder  radius,  and 


t  -  Gr/Rez 


(A3) 


In  equation  (A3),  Re  is  the  Reynolds  number  and  Gr  the  Grashof  number, 
given  by 


Re  -  u.  a/i/„ 

(A4) 

Gr  -  0ga3(Tw  -  T„)/v„2 

(A5) 

Hence,  buoyancy  effects  are  significant  for  large  overheats  and  small  flow 
speeds.  Applying  the  results  of  Yao  directly  to  the  present  geometry  with 
a  50°C  overheat  and  a  flow  speed  of  5  ms’1,  then  an  approximate  calculation 
gives 

a  -  0108m 
v.  -  108  x  10‘6  mV1 
p  -  5  X  10'*  "C-1 

Re  -  5  x  10s 

Gr  -  1-5  x  109 

t  -  6  x  10'3 
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The  above  calculation  suggests  that  buoyancy  effects  are  negligible  for 
x  <  0-14m  and  small  for  x  <  0-40m.  The  total  streamwise  extent  of  the 
heated  region  is  0-314m  and  furthermore  the  maximum  overheat  only  extends 
over  a  relatively  short  region.  Hence  buoyancy  effects  are  probably  at 
worst  very  small  in  the  present  experiments.  However,  if  it  was  required 
to  move  the  transition  location  much  further  back  in  a  future  experiment, 
then  buoyancy  effects  might  become  more  important. 
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PRESSURE  COEFFICIENT, 


HOT-FILM  GAUGES 
(AT  DIFFERENT  CIRCUM¬ 
FERENTIAL  LOCATIONS  ) 
FLUSH  WITH  SURFACE 


FIG.  1  TEST  MODEL  GEOMETRY  AND  PRESSURE  COEFFICIENT 
(FOUR-TO-ONE  ELLIPSOIDAL  NOSE) 
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(b)  WAVE-LIKE 
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FIG. 3  TYPICAL  HOT-FILM  SIGNALS  FOR  VARIOUS  FLOW  REGIMES 
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FIG.  5.  THEORETICAL  VARIATION  OF  SURFACE  TEMPERATURE  USING  DIFFERENT 

HEAT-FLUX  PROFILES 
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FIG.  13.  CORRELATION  OF  THEORETICAL  AND  EXPERIMENTAL 
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FIG.  14.  EFFECT  OF  SURFACE  HEATING  ON  TRANSITION 
LOCATION  (THEORY  AND  EXPERIMENT) 
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